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Abstract 

Models with an extra dimension generally contain background scalar fields in a non-trivial con- 
figuration, whose stability must be ensured. With gravity present, the extra dimension is warped 
by the scalars, and the spin-0 degrees of freedom in the metric mix with the scalar perturbations. 

o , 

CS| , Where possible, we formally solve the coupled Schrodinger equations for the zero modes of these 

g : ...... w.. . . .... . „ . . 

^ ■ superpotential, we are able to fully solve the system. We show how these zero modes can be used 

to construct a solution matrix, whose eigenvalues tell whether a normalisable zero mode exists, 
and how many negative mass modes exist. These facts are crucial in determining stability of the 
D ■ corresponding background configuration. We provide examples of the general analysis for domain- 

wall models of an infinite extra dimension and domain-wall soft-wall models. For 5D models with 
^ , two scalars constructed using a superpotential, we show that a normalisable zero mode survives, 

\^ ' even in the presence of warped gravity. Such models, which are widely used in the literature, are 

' therefore phenomenologically unacceptable. 
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I. INTRODUCTION 



A plausible way to extend the standard model is to embed it in one or more extra di- 
mensions. This opens up a new set of model-building tools which can help to solve a diverse 
range of theoretical and phenomenological problems, as well as yielding distinct collider sig- 
natures such as Kaluza-Klein (KK) modes. Almost all extra-dimensional models require one 
or more background scalar fields in some non-trivial configuration. For example, to generate 
a domain- wall which localises chiral fermions [l], to stabilise the size of a compact extra 
dimension , to generalise the Randall-Sundrum warped-space 1^ to a smoothed-out 
version jsl, or to cut off the extra dimension at a singularity jol, [ll. Domain- wall mod- 
els, whether they have an infinite or compact 9] extra dimension, make heavy use of 
background scalar configurations as a field-theoretic substitution for fundamental branes. 

Given the ubiquity and necessity of background scalar fields, it is important to under- 
stand both their statics and dynamics. The problem is best thought about in terms of a 
ground state, upon which exist perturbations. We want to know which particular scalar 
configurations have the lowest energy and are stable, and what the perturbations about a 
background lead to in terms of effective 4D modes. These two issues are closely related. 
The existence of negative-mass modes (tachyonic KK modes) signals an instability of the 
corresponding background. Massless modes may also signal an instability js], or be harm- 
less in the case of a translation mode. Positive-mass modes always exist and their precise 
spectrum is what distinguishes these extra-dimensional models at a particle collider. 

For the case of a single scalar field in a flat compact extra dimension, a general method 



for determining the lowest energy configuration has been worked out [10|,lll|- The inclusion 
of gravity in the analysis presents some complications because of the coupling of the scalar 
fields to gravity. This coupling generically warps the extra dimension p] and the scalar 
perturbations mix with the spin-0 degrees of freedom in the metric. For this warped case 
with one extra dimension there have been some general stability analyses with a single 
background scalar Q, [3], and some initial work on the multiple scalar case [o], [l^. For the 



case of a single scalar in multiple extra dimensions it has also been shown that the scalar and 
metric spin-0 modes mix [isl. Related analyses determining the spin-0 spectrum of multiple 
scalars in 5D have been done in the context of the AdS/CFT correspondence 16l-ll8l|. In 
particular, an algorithm for computing the scalar spectrum in general 5D compact models 
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has been prescribed [19|. 

Despite the comphcations introduced by gravity, it is still possible to find the effective 
coupled Schrodinger equations which describe the KK modes of multiple scalars in a warped 
background [9]. It is the aim of the current paper to, when possible, formally solve this set 
of Schrodinger equations for the massless KK modes — the zero modes — in the case of a 
5D bulk with no gravity (flat) and with gravity (warped). In addition to providing closed 
form solutions for the zero modes for a large number of cases, we shall also discuss how these 
solutions can be used to determine if any normalisable zero modes exist and whether or not 
the background is perturbatively stable. Some relevant examples shall be provided. 

One reason for studying the zero modes of the system is that, if they exist, they should 
play a large role at low energies in the effective 4D theory. For example, zero modes of a 
5D fermion are generally used to implement the fermions of the standard model 8j. When 
constructing domain-wall models using scalar background flelds in flat space, one always 
obtains a spin-0 zero mode corresponding to the broken translation symmetry. This degree 
of freedom is not welcome in the effective 4D theory since we have never observed such 
a particle, which would manifest as a "flfth force". Adding warped gravity can cure this 
problem since this removes the translation zero mode, as shown by Shaposhnikov et. al. 20| 
for the case of a single background scalar. In this paper this result is extended to the 
case of multiple background scalar flelds: the zero mode of translation does not survive, no 
matter how many scalars. It is possible though, and we shall give some explicit examples, 
that additional scalars introduce additional zero modes which do survive in the presence of 
gravity. Our examples of such models are constructed using the superpotential approach, 
and we argue that these models are phenomenologically unacceptable. 

The paper is organised as follows. In Section [TTl we flrst present the zero mode solutions, 
for the flat and warped cases, with general potential V and also specialising to a superpo- 
tential W with N = 2 scalars. For this latter case we give all four independent zero modes in 
closed analytic form. Section UTTl discusses the construction of a solution matrix, and how its 
eigenvalues can be used to flnd normalisable zero modes, and count the number of normal- 
isable negative modes. Following this we look in Sections IIVI and |V] at speciflc domain-wall 
models in both an inflnite and a compact extra dimension, and show that zero modes can 
survive in the presence of gravity. We conclude in Section IVI[ Appendix |X] summarises the 
method of reduction of order for ordinary differential equations. 



3 



II. THE ZERO MODE SOLUTIONS 



In this section we study spin-0 perturbations of real scalar fields in a 5D bulk for 
both a flat and warped extra dimension. The scalar fields are assumed to have some non- 
trivial background profile along the extra dimension, such as a kink. The aim is to derive 
formal solutions for the extra-dimensional profile of the zero modes, that is, the massless 
perturbations around the background. We shall concentrate mainly on the N = 2 case, and, 
for part of the analysis, specialise to potentials V that are generated by a superpotential W. 

Throughout this paper we work with the matter Lagrangian 



where qmn is the 5D metric with signature ( — + + ^i{x^, y) are scalar fields indexed 
by i = 1 . . . A^, is the 4D sub-spacetime, y is the coordinate of the extra dimension, 
and V{(^i) is the scalar potential. Repeated scalar field indices are always summed over. 
Perturbations of the scalar fields around some arbitrary background are written as 



The background profiles (t)i{y) depend only on the extra dimension, while the perturbations 
are general functions of all spacetime coordinates, and are required to be relatively small: 

iPi < 0i. 

From the equations of motion for the scalars one can obtain coupled, one dimensional, 
time-independent Schrodinger-like equations for the perturbations.^ The independent vari- 
able of the these equations is the extra dimension y and the eigenvalue corresponds to the 
mass of the KK mode of the perturbation. Where possible, we shall formally solve this 
system of coupled Schrodinger equations for the case of a zero eigenvalue. Since we have A^ 
second order, linear, ordinary differential equations (ODEs) we expect to obtain 2N linearly 
independent solutions. These 2N solutions are not necessarily physical (that is, are not 
normalisable) , but they do form a basis from which one can construct the unique solution 
for any given initial, boundary and/or normalisability conditions. The zero mode solutions 
are useful for studying the perturbative stability of the background configuration, as shall 
be discussed in Section UTTl 

^ In the case with gravity, we present + 1 Schrodinger equations plus one constraint equation, yielding 
effectively iV Schrodinger equations. 




(1) 




(2) 
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A. Flat case 



We analyse first the gravity-free, flat-space scenario where the action is given simply by 
S = f Cm d'^xdy. The discussion is divided into the case of a general V, and the case where 
V is generated by a (fake) superpotential W. 

1. General V 

In flat space, the general background equations are 

0'/ -V, = 0, (3) 
and the coupled Schrodinger-like equations for the perturbations are 

- <^'- + Vij<^j = n<^i . (4) 

Prime denotes derivative with respect to y, subscripts i and j on V denote a derivative 
with respect to $i and $j, and □ = d^dfj,. In addition, Vij, which is a function of 
must be evaluated on the background solution, $j = 4>i{y), where (pi{y) is a solution to 
equation Q. As usual, separation of variables of and y is the correct way to proceed 
here. For the sake of reducing the number of field variables we shall abuse notation slightly 
by using ipi to denote both the full perturbation which is a function of x'^ and y, as well as 
the separated factor that depends only on y. Separation of variables then proceeds as per 
ipi{x'^,y) = (fi{y)p{x'^), with p the 4D KK mode with mass m, such that Dp = m^p. 

One now solves equation (jl]) for the profiles ipi{y) and corresponding allowed mass val- 
ues, obtaining a tower of modes. Then one takes the original 5D action and substitutes 
$j = (f)i{y) + ipi{y)p{x'^) with (fi{y) a particular solution. The extra dimension can then be 
integrated out, leaving an effective 4D action for the mode p. To second order in p, this 
action is 



S= d'^x 



-e^., + Af(^-ld^pd,p-^m'p'^ 



+ (surface terms) . (5) 



Here, e^g is the energy density of the scalar background configuration. The normalisation 
constant for the mode is A/" = / ^pfdy and, so long as we pick a solution ipi that is normal- 
isable, one can scale said solution to obtain A/" = 1. The action then describes a canonical. 
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4D scalar field. The surface terms in the action are of the form J S'd'^x dy, where 5* is one 
of 

^1 = -(f)[ipi , (6a) 
-^2 = ■ (6b) 

The requirements that these terms independently vanish on the boundaries of the extra 
dimension, and that JV is finite, pick out the physical modes of the KK tower. 

As mentioned previously, we are interested in the zero modes, and shall look for formal 
solutions to equation (jlj) when Oipi = 0. For general V and all one solution is^ 

= 0: , (7) 

where 0j is the background scalar field configuration. Note that this solution is actually a 
vector of length N. It is the well-known translation mode of the background, since it is the 
first term in a Taylor expansion of the shifted background: (f)i{y + e) = (j)i{y)+e(j)i{y) + 0{e'^). 
For the case of = 1, where there are 2N = 2 independent solutions, we can use the y^-^^ 
solution to perform reduction of order (see Appendix [A]) and obtain the other solution: 

^r^'^=<P'l{<Pr'dy. (8) 

For A^ > 1 we could also perform reduction of order, but we would still have a relatively 
large (at least third order for N = 2) ODE to solve. At this point we shall be content with 
having only the translation solution for general A^. 



2. V generated by a superpotential 

It is possible to make progress and obtain an additional zero mode solution when the 
form of V is restricted to 

V = lw^, (9) 



^ Equation ^ is linear in (pi and we are free to scale any solution by an arbitrary constant, a constant 
which in some cases is dimensionful. For brevity, and because we are providing just formal solutions, we 
leave this constant out. Hence the units of this equation, and some of the equations that follow, do not 
match. 
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where W{^i) is a fake superpotential (it is just a formal construction and does not indicate 
a supersymmetry) . As it does on V, a subscript i on W denotes a derivative with respect 
to $i. In this case the background equations can be written as 

^[ = W,, (10) 

and the perturbation equation (jlj) becomes 

- + {W,,W,k + W,,uW,)vk = . (11) 

Here, Wij are to be evaluated on the background solution, and then they become functions 
of y. The utility of the superpotential approach comes from the fact that this perturbation 



equation can be factorised as [21 1 



{dySij + Wij){-dy5jk + Wjk)Vk = Oifi, . (12) 

Now when we look for zero modes, half of the solutions can be obtained by solving the much 
simpler equation 

i-dyS,, + W.,)cp, = . (13) 
For all N the translation mode still exists, 

^i'^ = , (14) 

and for = 1 the second solution will be given by equation ([8]). For N = 2 the situation 
becomes more interesting than the general V case. We can use the (pf'^ solution to reduce 
the order of equation f[T^ from two to one, and then solve the resulting first order ODE to 
obtain a second zero mode solution: 

^?' = W. I^dy + a.^. (15) 

Here we have defined 



J = exp 



{W^, + W22)dy 



(16a) 

X = a2Wi - aiW2 , (16b) 
Z = aia2{Wil - W2I) - {a\ - al)Wi2 . (16c) 

ai and 02 are constants and must be chosen such that the entity X (which is a function of 
y) is non-zero throughout the entire domain of y. For example, for systems where 0i is odd 



and 02 is even one can choose ai = 0, 02 = 1- We stress that different values of the do 
not generate independent zero mode solutions. 

At this point we need to make a few remarks about integration constants. There are two 
integrals in the solution (p] . The constant coming from the integral in J yields an overall 
normalisation factor for the zero mode solution. The constant in the integral in the first 
term in equation (IT^ pulls out a constant multiple of Lpf \ which effectively adds a multiple 
of this other zero mode solution. Thus our two integration constants amount to taking 
linear combinations of two linearly independent zero mode solutions. Alternatively, one can 
fix these constants of integration to zero and take linear combinations of equations ( IT^ 
and f lTSj) . Either way, we have a closed form for the general solution to equation f[T3|l when 
A^ = 2. 

There are two more linearly independent zero mode solutions for the N = 2 case. We 
cannot obtain them in closed form like the first two, but we can make some progress. Using 
the two known solutions ipf''^ and (ff^ the fourth order equation ( 1T2|) can be reduced to a 
second order ODE. This allows us to write the third and fourth zero mode solutions as 

V^f''^ = ^V^inA2^f , (17) 

where 

Ai = j J (wsLpf^ - W4V5f ^) dy , (18a) 

^ = y J {-W3ip^^ + w^if^^^^ dy , (18b) 

W3 = -Wuw-i - W12W4, , (18c) 

= -1^12^3 - W22W4 . (18d) 

The equations for w^^y) and w^^y) constitute the second order ODE which can be solved 
only with specific information about W. Since it is second order, there will be two sets of 
solutions for the pair {^3,^4}, and substituting these solutions into the equations for Ai 
and A2 will yield, through equation ( IT7|) . the final two independent zero modes. There are 
four integration constants in the above system of equations, as expected. Those in Ai and 
A2 add, respectively, a constant multiple of (ff^ and (ff^ to ff''^''. The other two come from 
solving for W3 and W4. (The constant from J can be absorbed in a rescahng of 1^3 and w^.) 
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B. Warped case 



We now repeat the previous gravity-free calculation for the case with gravity. It turns out 
that the Einstein constraint equation allows one to obtain additional zero mode solutions. 
The 5D action for N scalar fields coupled minimally to gravity is 

S = j d'xdy^(^-^R + C,n^ , (19) 

where = and M is the 5D Planck mass. Einstein's equations arising from this 

action are Gmn = '^i^'^Tmn where the stress energy tensor is Tmn = dM^d^^ + gMN^^m- 
We restrict our analysis to a warped metric ansatz, which is actually the most general 5D 
metric that respects 4D Poincare invariance, and is used extensively in realistic models. As 
for perturbations of the metric, we need only consider scalar perturbations, as vector and 
tensor perturbations decouple from the spin-0 sector 
the metric ansatz is 



With scalar perturbations F{x^, y), 



ds' = e-'"(l - 2F) T]^^ dxf'dx'' + (1 + 4F) dy' . (20) 

The warp factor exponent is (j(?/) and rj^i, is the 4D Minkowski metric. For consistency of 
small perturbations we require F ^ 1. The perturbations of the scalar fields are as in the 
previous section, equation ([2]). 

We now look for formal zero modes of this set-up, first in the case of a general scalar 
potential V, then in the case of V generated by a superpotential W. 



1. General V 



With a general potential V the background fields (piiy) and a{y) satisfy the equations 




(21a) 



(t)'l - Aa'cj)', -Vi = Q. (21b) 

By x^ scaling invariance, we are free to choose cr(0) = 0, leaving 2N integration constants 
for this set of equations. One of the redundant Einstein's equations, which is sometimes 
useful, is a" = n'^cpf. For the rest of this section, a and will be used to denote solutions 
to these background equations, and the potential V and its derivatives with respect to $i 
are to be evaluated on this background. 
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For the perturbation equations, it is best to work with the new variables xi^^iU) ^"^^ 
ipi^x^.y) defined by 

F = -^e'^x , (22a) 
= e^'^iJi . (22b) 

We shall write these + 1 components as a vector \1/ = {x, ipi) when it makes things neater; 

is an indexed version with m = . . . so that \l/o = x ^-^^d = ipi. In physical y 
coordinates, these perturbations obey the Einstein constraint equation (a detailed derivation 
can be found in Q|) 

X' -V2k<P[^, = 0, (23) 
as well as the coupled second order equation 

„ / 2a" 2x/2k(-(tV' +0'O \ o 

V v^, v^,; U„ = e2™^. (24) 

\2V2K{-a'(l)[ + C) (4a'2 - 2a")Sij + 6/cV>; + VijJ 

As in the fiat case, we again perform separation of variables, slightly abusing notation: 
F{x'^,y) = F{y)p{x^) and (pi{x'^,y) = (pi{y)p{x^). Equations fl23l) and f l2i|) allow us to 
solve for the KK tower of spin-0 modes, with extra-dimensional profiles F{y) and ifi{y)- 
Substituting these solutions in the metric and original field variables $j, computing the 5D 
action, and then integrating out the extra dimension yields the 4D effective action for the 
mode p: 

5 = y (fxN i^^d'^ pd^jtP — ^w?p^ + (surface terms) , (25) 
where the normalisation is 

A/- = / {X' + dy = j e-^- (^F^ + dy . (26) 

In deriving the effective action, we encounter three independent surface terms of the generic 
form f S'd^x dy, where S is one of 

^0 = Ae-^V , (27a 

1 



^1 = -^e^'^ie-'^'xy = -^e~'^{x' - Aa'x) , (27b) 
S2 = =^e'"'x{e~'"'xy - le-^^Ue'^i^iY = =^xix' - 28a'x) - + 2aV.) • (27c) 
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The subscripts here correspond to the order of perturbation. The last two equations in 
terms of physical variables are 

S, = ^e^^ie-'^'^Fy = ^^e-'^^iF' - 6a'F) , (28a) 



S2 = -^e^'^Fie-^'^^Fy - ^e-^>,^:, = -e-'^ (^^FF' - ^a'F^ + 1^,^ 

Sq must vanish on the y boundary for a background configuration to be physical. When 
looking for physical modes of perturbation, the solutions F{y) and ipi{y) must be such that 
A/" is finite and 6*1.2 vanish on the boundary. 

Let us now look for zero modes of this system, that is, when = 0. For this special 
case the constraint equation (1251) can be combined with the first row in equation to 
solve for x- 

X = ^ m + (2av: - 0- )^.] • (29) 

Thus the zero mode system is really just coupled, linear, second order ODEs. Ignoring 
finite normalisability and the vanishing of the boundary terms, such a system has 2N linearly 
independent solutions. Two of these solutions are 



= " , (30a) 



1 _^-2a 




V / 

where 

Bi= [ e-^'^dy . (31) 



These solutions were first derived by Shaposhnikov et. al. [20|] for the = 1 case (see their 
equation (3.6)), but the straightforward generalisation to all is also a solution. 

For = 1, \E'^^^ and \E'^^^ are the two linearly independent zero modes. For A^ > 1, we 
can use these known solutions to reduce the order of the system by two. We shall do this 
explicitly for the N = 2 case. Begin by using equation fl2^ to eliminate x the set of 
equations flMj) . This gives two second-order equations for tpi and 1^2- Now write this as four 
first-order equations and reduce the order by two using the method outlined in Appendix Rl 
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In terms of solutions, f{y), of this reduced ODE, the final two zero modes are 



/ 

where 













+ 




V ^^'2 ) 





G2 



G2 dy 



Hif+ / H^fdy 



i^i = 2/.V-p[log(0;e-)]', 
H2 = e^"^ {Q^^^ct>\ct>'^ + ^12) - 



'2ct J,/ 



The second order ODE that the auxiliary variable f{y) must solve is 

- /" + H^f + HJ = 0, 

where 



^3 = 

dy 



2k 



2 v^i^ 



0'i 



i74 = 4a'2 - 2a" + V22 - ^^12 + (^2a' - ^ 1 ^3 



(32) 



(33a) 
(33b) 
(33c) 
(33d) 

(34) 

(35a) 
(35b) 



The two independent zero mode solutions, \E'''3) and correspond to the two independent 
solutions for /. This system has four integration constants in total. Those in Gi and G2 add 
to \1>(3'4) a constant multiple of ^^^^ and ^^"^^ respectively. The other two constants come 
from the solution for /. 



2. Fake supergravity potential 



In the fake supergravity formalism 22, l23| the scalar potential is generated by a super- 
potential W: 

V = ^W^ - 2kV . (36) 
The background equations are then first order: 

a' = K^W , (37a) 
<^, = W,. (37b) 
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As for the general V case, we work with the variable \& = {x, ipi) for the perturbations. The 
Einstein constraint equation is 

x'-v^KVTi^i = 0, (38) 
and the equivalent of equation f l24|) factorises to give 



{dy + U){-dy + U)^ = e^''U^> , (39) 

where 

U=\ \ (40) 




We now look for zero mode solutions to equations f l38|l and fl39|) . To begin with, we have 
the two solutions found in the general V case for all A^, written here in terms of W: 

= r , (41a) 



(41b) 



We now concentrate on the N = 2 case and perform reduction of order on the system of 
equations. Due to the factorisability of the perturbation equation fl39|l . we proceed here in a 
different manner than we did in the case for general V. We begin with the third order system 
{—dy + f/)^ = and use the two known solutions \['(^'^) to reduce the system to a single first 
order ODE, which we solve for the third solution \E''^^\ To get the fourth solution, we take 
the full sixth order system {dy + U){—dy — f/)\E' = and eliminate x using equation fl2^ 
(using W instead of the background fields). We then use solutions vl/^^'^'^) to reduce the 
resulting system from order four to order one. This final ODE can be solved to find '^^^K 
The two additional solutions are 

^(3) ^ ^^^(1) ^ ^^^(2) ^1 I ^ (42a) 

^ re ^"^a,- 



1 ^2cr W2 
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The auxiliary factors are 

C2 
Di 
D2 
D3 
T 



V2J 



-2a 



X 
JZ_ 

BiT 



WZ\ 



BiC[ ) dy 
1 



J 



dy, 



D',i-Ci+T)dy 
W 



Wo 



-{W[ - K^WWi) 



1 Wu 2 W 
h K ttj 

2Wi X 



(43a) 
(43b) 
(43c) 
(43d) 
(43e) 
(43f) 



W'Wi 

Note that everything here is ultimately defined only in terms of W , its derivatives with 
respect to $, and a, all evaluated on the background. Once W is given, everything else can 
be computed in a closed form, including the four linearly independent zero modes (for the 

= 2 case). Also note the identities W = Wf + and W[ = W^Wi + W12W2. 

It is not immediately obvious, but there are only three independent integration constants 
in the definition of and four in These constants pull out constant multiples of 

lower zero mode solutions. In effect, '^'^^^ is the most general zero mode solution. 

There are two conditions that allowed us to find the general zero mode solution in closed 
form for = 2 in the fake supergravity case. One, there is a constraint equation, and two, 
the rest of the perturbation equations factorised. In contrast, for the fiat case with W we 
did not have the constraint equation, and for the warped case with general V we could not 
factorise. 

The full zero mode solution that we derived for general V in the warped case, equa- 
tion (1321) . is equivalent to the solution '^'^^^ found in this section, although they are written 
in manifestly different ways. It is straight forward to write them in equivalent ways, allowing 
us to find the solutions to equation flMl) for /: 



yd) 



Wi ^ 



3 ; 



. w 

-^"——dv 



(44a) 
(44b) 



(The first of these was found by inspection, the second, by reduction of order using the 
first.) Putting these solutions in equation (1521) . along with the relevant substitutions for 
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the backgrounds a and 0j in terms of W ^ yields equivalent expressions for the zero mode 
solutions \Ef(^'^). These can be used in place of equations (142 ap and fl42bp if desired. 

Unfortunately we cannot use these solutions for / to intelligently deduce the correct 
solutions in the general V case. This is because J appears in /, which is computed from 
W\\ and 1^22- These latter functions cannot be written in terms of V ^ its derivatives, and/or 
the fields a and 0j. We also remark that while one can recover the known fiat case zero 
mode solutions by taking /t — > in the warped solutions this does not produce any new 
solutions for the flat case. 



C. Summary of zero mode solutions 

Let us recall the main results of this section. For scalar fields in a flat and warped 
extra dimension there exist IN linearly independent zero mode solutions for perturbations 
around a background configuration. These formal solutions may or may not be physical; 
physicallity is obtained by demanding the normalisation is finite and the surface terms 
Si vanish at the boundaries of the extra dimension. 

For = 1 scalar field, the two zero modes are completely determined for general V ^ for 
both the flat case, equations ([7]) and ([8]), and warped case, equations fl30ap and fl30bp . 

For N = 2 scalars we have derived the following results: 

• Flat space, general V: one explicit solution, equation ([7]). 

• Flat space, using W: two explicit solutions, equations f lT^ and f lTSj) . and a second 
order ODE for the full solution, equation (fT7|l . 



Warped space, general V\ two explicit solutions, equations f l30ap and fISObp . and a 
second order ODE for the full solution, equation ( 132|) . 



• Warped space, using W: four exphcit solutions, equations (14 lap . (I41bp . (142 ap 
and fl42bll . 

For > 2 we can say the following for both general V and a V generated by a super- 
potential W . In a flat extra dimension the zero mode of translation always exists, given 
by equation ([7]). Extra dimensions that are warped admit two closed form solutions, equa- 
tions (ISUip and (!30b|) . 
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III. THE USE OF ZERO MODES 



By knowing the formal zero mode solutions of a system (they need not be physically 
normalisable) , we can deduce some important physical properties of that system. We shall 
provide some practical remarks on using zero modes solutions to 

• look for linear combinations that give normalisable, physical zero modes; 

• check for perturbative stability. 

Our discussion here is restricted to systems whose background configuration has definite 
parity, that is, (piiy) are either even or odd under y —y (different fields can have different 
parities).^. For the warped case, a must always be even, due to the Einstein equation 
a" = K,'^4>f. Furthermore, since we choose a{y = 0) = and also demand that at least one 
of the scalars have (f)'^{y = 0) 7^ 0, we have that a is strictly monotonically increasing. Then, 
without negative tension branes, the extra dimension can only end when cr — ?► 00 (otherwise 
the junction conditions coming from Einstein's equations cannot be satisfied at the patching 
points). We then have two scenarios: either a diverges only as ?/ — t- 00 and we obtain an 
infinite extra dimension j^, or a and at least one scalar diverge at some finite y value and 



we get a soft wall js, 7|. Examples of both of these types of spaces will be presented in the 
following sections. 

The 2N linearly independent zero modes of a system can be written in an infinite number 
of ways, as they form a basis for the set of all solutions to the massless perturbation equation. 
Regardless of the linearly independent solutions that are obtained, one should be able to 
compute characteristic, physical properties of the system in an unambiguous way. To find 
these characteristics, our idea is to construct a specific matrix of the zero modes solutions 
(which will be a function of y) and then compute the y dependent eigenvalues of this matrix. 
Looking at these eig envalues is an extension of the idea of looking at the determinant of the 
solution matrix 
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2^. 



Specifically we want to construct an x square matrix whose columns are vectors 
of formal zero mode solutions, that is, a single column is a vector whose N entries are a 



^ Note that even though the backgrounds (pi must have parity, there is no restriction on the fuh field $i 
nor the perturbations. 
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particular solution (pi^ Actually, we want to construct two such matrices: a matrix ME{y) 
for the even solutions — those that have </7j the same parity as the corresponding background 
field (f)i — and Mo{y) for the odd solutions, whose parity is opposite the background field. 
The initial conditions (values of the perturbation Lpi at y = 0) of the solutions in the Me 
(Mo) matrix will form a basis for an arbitrary even (odd) mode. For N scalars there will 
be linearly independent even and odd solutions, so our matrices will have N columns. 

As long as the above criteria are satisfied, the actual initial conditions of the Me,o 
matrices are not important. But for clarity we shall describe a simple realisation. Let the 
system have A^ background fields with parity Pi G {0, 1} such that an even (odd) field 
has value 1 (0). Then the even solution matrix has initial conditions 



ME{y = 0) 



Pi 

P2 



\ 



yO ... Pnj 



( 



M'Eiy = 0) 



i-Pi 
I-P2 














1-Pn 



(45) 



y u u . . . ± — j^jv J 

The odd solution matrix has opposite initial conditions: Mo{y = 0) = M'^{y = 0) and 
M'ply = 0) = ME{y = 0). Given these initial values, one must then compute all entries in 
the two matrices as a function of y, up to the boundary of the extra dimension. This can 
be accomplished by using the closed form expressions for the zero modes in the previous 
section, or by directly integrating the coupled ODEs describing the perturbations. In the 
former case the integration constants in the closed form solutions must be chosen to achieve 
the correct initial conditions. In the latter case the initial conditions in Me,o can be used 
directly as initial conditions in, for example, a numerical ODE solver. 

The full ME,o{y) matrices form a basis of even and odd solutions because their initial 
conditions form a basis of all possible initial conditions. The set of all zero mode solutions 
is generated by the matrix products Me ■ ol and Mq ■ ol, where a — (ai, . . . , ctjv), ctj e IR is 
a vector of constant coefficients. 

We can now compute the eigenvalues of Me,o- They will be functions of y and there 
will be A^ of them; let us denote them by \f''^{y). These eigenvalue functions give us a 



* Recall that for the warped case the gravitational perturbation x could be solved for in terms of il^i, so only 
the degrees of freedom, equivalently the <^i, are needed to construct our matrix. Thus, the following 
discussion is valid for both the flat and warped scenarios. 
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lot of information about the stability of our background configuration. Assuming that the 
condition for normalisability is that a perturbation must vanish at the boundaries of the 
extra dimension, at y*, we make the following two conjectures: 

1. For each Xf ^ as y ^ there exists a corresponding, normalisahle, even zero 
mode given by Me ■ oti, where ai is the eigenvector for \f{y*)- The correspondence is 
one-to-one: the existence of a normalisahle mode implies the vanishing of one of the 
eigenvalues at y^. An equivalent statement is true for the odd sector. 

2. The number of negative mass modes in the full spectrum of perturbations equals the 
number of times the eigenvalues ' pass through zero in the domain < ?/ < . 

We shall sketch the proof for the first conjecture. Let M be either Me or Mq- A 
normalisable zero mode exists if we can find a linear combination of formal mode solutions 
that vanishes as y — y*.^ That is, we want to find a constant, non-trivial vector a such 
that M ■ a — )■ for ?/ — i- This is simply an eigenvalue equation with a zero eigenvalue, 
and with eigenvector a. Thus, we want to find the eigenvalues of the solution matrix 
M, called Xi{y), and we want at least one of these eigenvalue functions to tend to zero at 
the boundary of y. If there exists such a Aj, then the corresponding eigenvector evaluated 
at ?/* gives the coefficients needed to construct a normalisable zero mode. Conversely, if a 
normalisable zero mode is known to exist, then one can find the vector a such that M -a = 
at ?/ = y*, and so M has an eigenvalue function which vanishes at the boundary. 

For the second conjecture we do not provide a proof. It is based on a closely related 
theorem given by Amann and Quittner [2J]. They work with a system of coupled radial 
Schrodinger equations and the wavefunction values must all vanish at the origin; effectively 
they are looking only for solutions where all wavefunctions are odd. Their proof should be 
adaptable to our second conjecture stated above, including correct handling of the weight 
function e^'^ in the warped case in equation (12^ . We do not attempt to construct the 
proof here, and the analysis in the following sections is largely independent of it. Our 
interest in presenting the second conjecture is to show that the eigenvalues Xf'^ contain more 
information than just whether or not a normalisable zero mode exists. For an application 
of Amann and Quittner's theorem to a system with 16 components see 25|. 



^ Since we are dealing with a Schrodinger-like equation, solutions at infinity either oscillate or behave expo- 
nentially. If a solution asymptotes to zero then it must decay exponentially, implying square integrability. 
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To summarise, we claim that the eigenvalue functions Xf''^{y) of the general solution 
matrices ME^oiv) give all the information about the perturbative stability of a specific 
background configuration, for both flat and warped extra dimensions. They tell the number 
of unstable modes, if any, and whether or not the configuration is critically stable, that is, 
it admits a normalisable zero mode. For phenomenological reasons, one generally tries to 
construct models that are free of zero modes. 

In the following sections we shall apply our first conjecture to some example models to 
show the existence, or lack thereof, of zero modes. 

Before moving on to the examples, let us discuss the massless translation mode associated 
with a background. For the fiat case, there is always a formal solution corresponding to 
translations of the background, equation ([7j). Assume this solution is normalisable. We 
would like to know what happens to this mode when one adds gravity to the system. For 

= 1 scalar, Shaposhnikov et al [2^ have shown that when gravity is turned on this mode 
is no longer massless and becomes instead a wide resonance. 

Now consider the existence of a translation mode for general with gravity. The solution 
vl/(i) given by equation f l30ap looks like it has the right form for such a mode, as the tpi 
components are exactly Even though the actual perturbations are ipi = e'^^ipi = e^'^^'j, 
this solution still has the correct initial conditions for it to be a valid translation mode: 
V9j(0) = 0'j(O) and fi{0) = 4>i{0) since a{0) = cr'{0) = 0. These initial conditions are enough 
to uniquely specify the mode solution, so '^^^^ is the solution with the initial conditions 
of a translation mode. But this mode is non- normalisable, since cr — )■ 00 as ?/ — ?■ y*.^ We 
therefore conclude that for general A^ with a warped extra dimension the usual translation 
mode is rendered non-normalisable by the introduction of gravity. 

From a phenomenological point of view the absence of a translation mode is good news, 
since massless spin-0 particles are not seen in nature. But, for A^ > 1, it may be that 
there are additional zero modes in the spectrum that survive when gravity is turned on. 
For example, there may exist a zero mode that both translates and dilates the background, 
and so has different initial conditions to '^^^^ and is normalisable. We show in the following 
sections that such modes do in general exist, and the addition of gravity does not in general 
remove all zero modes from the spin-0 particle spectrum. 

^ Unless there exists such that e^'^(f>'^ — ^ as y — > y, , with cr" = n^cj)'^ ■ 
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IV. DOMAIN- WALLS IN AN INFINITE EXTRA DIMENSION 



Our first example model has an infinite extra dimension with N = 2 scalars in a kink- 



lump configuration. This type of set-up was used in js] in an attempt to realise the standard 
model confined to a domain-wall. We shall present two incarnations of this model. The first 
is effectively that presented in 8| and is described by a straightforward quartic potential; 
we shall call it the \^-model. As will be shown, in the fiat case this l^-model has a single 
zero mode, whereas in the warped case the zero mode is absent. The second incarnation is 
based on the fake supergravity approach and is described by a superpotential W; we call this 
model the PF-model. We shall see that although the background solutions are qualitatively 
the same as those in the V^-model (and in fact can be made exactly the same for certain 
choices of parameters) the W^-model posses two zero modes in the fiat case, one of which 
survives when gravity is turned on. 

A. Kink-lump model in flat space 

The Lagrangian of the fiat-space l/-model is given by equation ([T]) with N = 2 scalars, 
$1^2, and potential 

2 2 4 

It is not our aim here to perform a complete analysis of this model. Instead, we shall use it 
to give a constructive example of how one looks for normalisable zero modes, and show the 
effect of adding gravity. 

We first restrict ourselves to the region of parameter space where all five parameters are 
positive, cf 2 — /i^ > and AlXv'^ — /it- > 0. This ensures that the kink-lump configuration 
is stable (has no negative modes) |26|. Next, we impose the relation 2(A — c)/U^ = (2cA — 
4/A — c^)v'^ which allows us to obtain analytic solutions for the background 8|: 

0i = ftanh(%), (j)2 = Acoshr^{ky) ^ (47) 

where fc^ = cv"^ — jj?^ and A? = (2yU^ — cv'^)/X. These are solutions of equation ([3]). The 
background configuration is straightforward with (pi the kink and 02 the lump; we do not 
plot these functions for the flat case, but see Figure [3] for qualitatively similar plots of 0i,2 
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FIG. 1: Eigenvalues of the zero- mode solution matrices ME^oiu) (left) and the normalisable zero 
mode of translation (right) for the flat-space kink-lump l/-model. The single odd eigenvalue that 
asymptotes to zero at large y signals the existence of the translation zero mode. Since the other 
three eigenvalues diverge, there are no other normalisable zero modes. 

in the warped case. Note that in all our plots we show only the y > half of the fields. The 
other halves are found by relevant parity transformations. 

Using this background, we now solve equation (jlj) for the four independent zero modes, 
two which are even and two which are odd, and construct the zero-mode solution matrices 
ME,o{y)- Recall that even and odd are relative to the background configuration, so, for 
example, the even set of zero modes in M^(y) have ipi odd and even. Parameters we 
choose are / = 0.7, v = 1.0, c = 1.5 and A = 0.4, with derived parameter /x^ = 0.98636. 
This choice gives typical looking solutions which exhibit the behaviour we are interested 
in. It is not a fine-tuned choice and small variations in the parameters give qualitatively 
similar results. Having computed the two 2x2 matrices ME^oiv) we then compute their 
two eigenvalues (four eigenvalues in total), which are plotted as a function of y on the left 
in Figure [U In order to accommodate the large range of the eigenvalues we implement 
a quasi- log scale by plotting arcsinh[Af ''^(y)]. 

This plot gives a quantitative summary of the stability behaviour of the background con- 
figuration of the V^-model. According to conjecture one, the existence of an odd eigenvalue 
A*^ with A*^ — 7- as ?/ — 7- oo, as is apparent in the plot, implies the existence of a nor- 
malisable zero mode. The eigenvector corresponding to A*^ evaluated at large y gives the 
linear combination of formal zero mode solutions which yield the normalisable zero mode. 



21 





.2 I I I I I I I I I I I I I I I I I 

012345678 012345678 
extra dimension y extra dimension y 



FIG. 2: Eigenvalues of the solution matrices (left) and the normalisable even zero mode (right) for 
the flat-space I^-model. There are two normalisable zero modes for this model, the other one is 
the odd translation mode and is exactly the same as in the l^-model. 

The resulting initial conditions for this normalisable mode are ((y9i(0), ^92(0)) = (1,0) and 
{(p\{^), v?2(0)) = (0) —0.77912), and the mode is plotted on the right in Figured] This mode 
is exactly the translation zero mode given by equation ([7]), as expected. The fact that the 
other three eigenvalues diverge at large y implies, by conjecture one, that there are no more 
normalisable zero modes for this particular background configuration. By conjecture two, 
since none of the eigenvalues cross zero there are no negative mass modes, which is again as 
expected due to our choice of parameters. 

Let us now consider the flat ly-model, which admits similar kink-lump configurations as 
the \^-model, but has different behaviour when it comes to the zero modes. The W^-model 
has a potential described by equation (|9]) with superpotential given by 

W = a^i-h^\-c^i^l. (48) 

This model always has the analytic kink- lump background solution given by equation P7|) . 
with the constants in the solution given in terms of the parameters in as f ^ = a/36, 
k'^ = Aac^/Sb and = a/c — 2a/3b. For the following analysis we choose parameters to give 
the same v, k and A as in the ^-model, namely a = 1.14017, b = 0.38006 and c = 0.35834. 
The outcome of our analysis is not so dependent on parameter choice since we simply want 
to demonstrating the existence of zero modes. 

As before, we compute the solution matrices and their eigenvalues, the latter of which are 
shown on the left in Figure [2J From this plot we see that the W^- model has two normalisable 
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zero modes, one odd and one even, and no negative modes. The normalisable odd mode is 
the translation mode, and is the same as in the \^-modeL The normahsable even mode has 
initial conditions (y9i(0), ^92(0)) = (0, 1) and {ip[{0) , ip'2{0)) = (—0.77912,0), and is shown on 
the right in Figure [2l Physically, this even mode expands both y and 02, or shrinks both. 
Although any linear combination of the odd and even normalisable zero modes is again a 
normalisable zero mode, what is important is that there are a total of two massless physical 
degrees of freedom. The choice between an odd and even basis, or some other basis without 
parity, will depend on the physical problem at hand. 

Note that in |^ it is shown that for models generated by a superpotential W, exciting the 
zero modes corresponds to changing the integration constants in the first order equations of 
motion for the background, (j)[ = Wi. In our example here we have two fields, two integration 
constants, and so two zero modes. The odd mode changes the value 0i(O) and the even mode 
changes 02 (0). These modes are massless because changes in the initial values 0i(O) do not 
change the energy density of the configuration. 

In summary, even though the \^-model and W^-model have the same background config- 
uration, the former admits only one normalisable zero mode, while the latter admits two. 
Since the VT-model contains some extra symmetries owing to its supersymmetric nature, it 
has the extra zero mode. These results are instructive, although not particularly profound 
in their own right. The reason we have constructed these two models is so we can compare, 
qualitatively, what happens when gravity is turned on and the extra dimension is warped. 



B. Kink-lump model in warped space 

We now look at the existence of zero modes for the kink-lump model when gravit y is 
turned on and the infinite extra dimension is warped, as per a smoothed-out version jsl. l2^ 
of Randall- Sundrum [J]. We shall analyse both the V- and VT-models. 

The action is given by equation f|T9l) with metric ansatz fl20|) . The scalar potential is 

K.arp = V^ + A, (49) 

where V is the flat space potential, equation fj46|) . and A is the bulk cosmological constant 
required to fine-tune flat 4D slices in the warped space. For stable solutions, constraints on 
the parameters in the potential are the same as for the flat space case. We can again obtain 
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FIG. 3: Background configuration for the kink-lump model with a warped extra dimension, 
analytic background solutions, but the relations are now different, being 
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where k"^ 



[cv 



;U^)/(1 + Ak^v"^). Recall that k is the 5D Newton's constant. With these 



relations the background configuration is 

a = K^v^ log [cosh(%)] , 01 = viaxih.{ky) 



^2 = V cosh ^{ky) . 



(51) 



A representative choice of parameters is k = 1.0, / = 0.7, v = 1.0 and c = 1.5, with derived 
parameters A = 0.2, /i^ = 0.41667. Plots of the background fields are given in Figure [31 In 
the gravity-free limit k — )■ 0, and we have a — t- while 01 2 retain their kink-lump structure, 
so we obtain similar solutions as in the fiat space case. But we must make it clear that the 
point in parameter space we have chosen for the warped case is not the same (but it is close 
to) the point in parameter space that we analysed in the fiat case. Nevertheless, since we 
are interested only in qualitative features of the set-ups, we can still make a fair comparison 
between the fiat and warped configurations. 

Given this warped background we can proceed to compute the zero mode solution matrices 
ME,o{y)- The metric perturbation x (and its counterpart F) can be written in terms of the 
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FIG. 4: Eigenvalues for odd and even zero modes for the model with a warped extra dimension. 
All eigenvalues diverge for large y so there are no normalisable zero modes. Since the eigenvalues 
do not pass through zero, so there are also no negative modes. 

ipi (counterparts y^j), so we only need the latter to construct the solution matrices. Now, if 
we use ipi to construct the matrices and look for eigenvalues that tend to zero for large y, 
we will obtain modes that are normalisable with respect to the integral j ip'^dy. But what 
we really want are modes that are normalisable as per equation fl2^ . So in fact we should 
construct solution matrices from y^j, which, by equation ( ]22b[) . is simply e^" ME,o{.y), where 
Me,o here is constructed from t/^j. 

Figure m shows the eigenvalues of the solutions for the warped V^-model. All of the eigen- 
values diverge for large y so there are no zero modes in the physical spectrum. In particular, 
the translation zero mode does not survive in the presence of gravity, in concordance with 



the results of 



20|. Similarly, since none of the eigenvalues pass through zero, there are no 



negative modes in the spectrum either. This is also as expected since we restricted our 
parameters so the configuration would be stable. In summary, the V^-model in warped space 
has no normalisable negative modes and no normalisable zero modes. 

Now consider the W^-model in warped space. We use the same superpotential as given 
by equation fHS]) . but now the derived potential V is modified, as per equation fl5Bl) . Such 
a potential is again, as in the fiat space case, qualitatively similar to the y-model. In fact, 
for the choice of parameters 6 = c, we can get analytic solutions of exactly the same form 
as the warped model, equation flSTj) (such a model is used in [27|). The parameters of 
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FIG. 5: Eigenvalues of the solution matrix (left) and the single normalisable zero mode (right) for 
the warped-space M^-model. Compare with the flat-space T^-model, Figure [21 which has two zero 
modes, and the warped-space y- model, Figured! which has no zero modes. 

this solution in terms of the parameters in W are = a /3b and k'^ = Aab/S. To obtain a 
background with exactly the same form as the one we used in the analysis of the V^-model 
we choose a = 0.69821 and b = c = 0.23274; the background is show in Figure [31 The 
conclusions that we shall draw regarding zero modes are generically the same for a large 
parameter range, but we make this choice so we can compare with the V^-model. 

The eigenvalues of the solution matrix for the physical perturbations for the warped 
ly-model are shown on the left in Figure [21 As can be seen, there is a surviving even zero 
mode. The initial conditions for this mode are (y9i(0), ^92(0)) = (0, 1) and {lp[{0) , {p'2{0)) = 
(—0.172148,0) and the mode is plotted on the right in Figure [51 It is normalisable, as per 
equation (1261) . and the surface terms, equation (128|) . vanish at large y. The mode is therefore 
present in the physical spectrum. It has a qualitatively similar form to the even zero mode 
in the fiat 1^-model. Finally, no eigenvalue crosses zero, so there are no negative modes in 
the spectrum, a result which is already known for general W j^. 

Even though the potentials of the four models we have looked at (fiat and warped, V- and 
H^-models) are qualitatively very similar and admit the same background configurations, 
they show very different behaviour when it comes to having zero modes in the physical 
spectrum. Our general conclusion is that adding gravity in the form of a warped extra 
dimension will remove the translation zero mode from the spectrum, but will not necessarily 
remove other zero modes. In this section we have also explicitly shown how the eigenvalues of 
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the solution matrices ME^oiv) allow one to easily find normalisable zero modes. Furthermore, 
our analysis of the warped ^-model shows, at least for some values of the parameters, that 
it contains no zero modes. It is therefore phenomenologically acceptable to use this type of 
set-up for constructing realistic models, as is done in j8|. 



V. DOMAIN- WALL SOFT- WALL MODELS 



In this section we briefly analyse another example model, one with a compact extra 
dimension and N = 2 scalars. The potential is generated by a superpotential with gravity 
and we show that a zero mode again survives in this compact set-up. The model was first 
presented in 9I] as a realisation of a domain-wall soft-wall model, where the extra dimension 
is dynamically compactified by the formation of curvature singularities. 

The superpotential is 



W = a sinh(i/ $1) + (a $2 - 6 $2 



(52) 



For background configurations of definite parity with both scalars odd there is a unique 
solution to the first order equations of motion: 



a 



log [cos (au'^y)] — |l + 41og 



cosh I V Sab y 



— arctanh 



tan 




tanh I V Sab y 



cosh 2 (^VSaby^ | , (53a) 

(53b) 

(53c) 



The edge of extra dimension is fixed at = njlav^. We choose parameters a = 1.0, 
u = 1.4, a = 0.5 and b = 0.3 and plot the background in Figure El In it was shown 
that enforcing odd parity on the fields $1,2 themselves, as opposed to just the background 
configuration, ensures that there are no normalisable zero modes in the spectrum. We shall 
now show that relaxing the parity condition leads to the appearance of a zero mode that 
de-stabilises the background. 

For the background configuration and choice of parameters given above we compute the 
solution matrices ME,o{y) for the physical perturbations ipi. The eigenvalues of these two 
matrices are show on the left in Figure [3 Three of the eigenvalues diverge as y — )■ y*, 
but the other one remains finite. Even though our first conjecture states that we must 
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FIG. 6: Domain-wall soft-wall background configuration, (pi acts like a dilaton, while (j)2 takes on 
a (very wide) domain-wall solution. The extra dimension ends at a physical curvature singularity 
at = 0.80143. 

look for eigenvalues that tend to zero to find normalisable modes, we find that this finite 
eigenvalue does actually correspond to a properly normalisable mode. Since our space is 
finite, if yjj and F remain finite throughout the extra dimension they will be physically 
allowed perturbations (a small multiple of them will be small compared with 0, and a). The 
initial conditions for the normalisable mode corresponding to the finite odd eigenvalue are 
(<^i(0),<^2(0)) = (-0.19847, 1) and (<^'i(0), (^2(0)) = (0,0). This mode is shown on the right 
in Figure [71 It is normalisable as per equation fl26|) . and the relevant surface terms in the 
effective 4D action vanish because e~^" — > at the boundaries of the extra dimension. This 
zero mode has opposite parity to the background, so if one does not enforce parity on the 
fields themselves this mode will exist in the 4D spectrum and the background configuration 
will not be stable. 

As we have shown, the eigenvalues of the solution matrices are also useful for finding zero 
modes when the extra dimension is compact and the perturbations are allowed to be finite 

ri 

over all y. The domain-wall soft-wall model in [9| contains a normalisable zero mode which 
must be removed by enforcing parity on the fields themselves (the model-building philosophy 
of these set-ups forbids adding fundamental branes to restrict the boundary conditions of 
the KK modes, thereby eliminating any zero modes). Alternatively, we suggest that it may 
be possible to construct a domain-wall soft-wall model using a normal potential V that does 
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FIG. 7: Eigenvalues of the solution matrices (left) and the normalisable odd zero mode (right) for 
the domain-wall soft-wall model. 

not have the extra symmetries inherent in the superpotential approach, and hence does not 
have a surviving zero mode. 

VI. CONCLUSIONS 

In this paper we looked at scalar perturbations around a background configuration, and 
were concerned with finding zero modes. Both flat and warped extra dimensions were stud- 
ied, with both a general scalar potential V , and a potential generated by a superpotential W . 
The results can be split into three main parts: analytic zero mode solutions, use of these 
solutions, and examples. In Section [IT] we presented analytic expressions for formal zero 
mode solutions to the perturbation equations, with particular emphasis on the N = 2 scalar 
case. For N = 2 scalars in warped space with a potential generated by a superpotential we 
found analytic closed-form expressions for the four, linearly independent zero modes solu- 
tions. Section [TTll discusses the use of formal zero modes, and here we made two conjectures 
which use the eigenvalues of the solution matrices: one states how to find normalisable zero 
modes, the other tells how many normalisable negative modes there are in the spectrum. In 
Sections IIVI and IVl we looked at some specific example models to demonstrate the workings 
of these conjectures, and to show that normalisable zero modes can survive when the extra 
dimension is warped. 

Our general conclusions regarding extra-dimensional model building are the following. If 
one uses a general potential V then in fiat space there will exist the translation zero mode. 
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which is removed from the 4D KK spectrum when the extra dimension is warped (for = 1 

n 

scalar this conclusion was made in [20|). But it is not always the case that all spin-0 zero 
modes are removed by the inclusion of gravity. We have shown that models which have N = 2 
scalars and generate the potential from a superpotential generally admit an extra zero mode 
which survives in the presence of gravity. Such superpotential models are widely studied 
in the literature for the reason that they give first-order equations of motion. But without 
some extra input, like a fundamental brane or forced parity, these superpotential models 
will be phenomenologically unacceptable due to the presence of massless spin-0 degrees of 
freedom, something which we have not observed in nature. 

Finally, the general solutions we have found for the fake supergravity scenario with N = 2 
scalars rely only on the metric ansatz (l20l) and should have wider applicability than to the 
domain-wall models that we emphasise here. For example, the inclusion of fundamental 
brane terms would change only the boundary conditions; the bulk solutions we have found 
will remain the same. Extending the zero mode solutions to more than one extra dimension 
may also be possible, following the analysis of jl5|. 
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Appendix A: Reduction of order of a set of linear homogeneous ordinary differential 
equations 

Any set of linear homogeneous ODEs which is n*"^ order can be easily recast as a set of 
n first order differential equations. We assume this has been done and that the resulting n 
dependent variables fi{y), where i = 1 . . .n and y is the independent variable, satisfy the 
equations 

m = A,{y)f){y). (Al) 



^ For related results on the survival of light (but not exactly massless) spin-0 states, see 
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There is an implicit sum over j. This equation has n independent solutions. If we know m 
of these solutions then we can perform reduction of order on f lAip to obtain a set of n — m 



coupled first order ODEs. We outline this procedure, which follows closely that given in 29 |. 

Write equation ( ]A1|) as a matrix equation, f'{y) = A{y) ■ f{y), where f{y) is a vector 
of length n made of the dependent variables. Let Fg be the n x n solution matrix of this 
ODE, such that each column of Fs is an independent solution vector f^'^\y). It must be 
that det Fg for all y. Then Fg is known as the fundamental matrix of the set of ODEs 
specified by A, since Fg determines A uniquely hj A = F!,F~^ (the converse is not true since 
Fg ■ (const matrix) is also a fundamental matrix of A). 

Reduction of order then proceeds as follows. Assume we know m columns of Fg (that is, 
m linearly independent solutions of the ODE) which we label f^'^^^ with i = 1 . . .m. Then 
construct the n x n matrix 

U = {f^''\ . . . , a^^+^\ . . . , a(")) , (A2) 

where the a'^^\ j = m + 1 . . .n, are linearly independent constant vectors. They can be freely 
chosen, so long as det U 7^ for all y. Usually the a^^'^ can be unit vectors. Now change 
variables to g by the definition f = U ■ g and the equation f = A ■ f becomes 

^?' = f/-i-A-(0,...,0,a('"+i),...,a("))-^?. (A3) 

The first m components of the vector g do not appear on the right-hand-side of this ODE 
(after multiplying the matrices out), so this procedure decouples m of the equations. Call 
the solutions to equation ( IA3I) g^^^\ We know m of these are just constant vectors: g^^^^ = 
(1, 0, . . . , 0), g^^'^^ = (0, 1, 0, ... , 0) and so on to g'^'^"^\ The remaining n~ m solutions are 
to be determined using other techniques. Once they are found, the remaining solutions to 
the original ODE are given by Z*^*^-* = U ■ g^^^\ 
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